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ABSTRACT 

We measure the temperature of warm gas at planet-forming radii in the disk 
around the classical T Tauri star (CTTS) TW Hya by modelling the H 2 fluores- 
cence observed in HST /STIS and FUSE spectra. Strong Lya emission irradiates 
a warm disk surface within 2 AU of the central star and pumps certain excited 
levels of H 2 . We simulate a ID plane-parallel atmosphere to estimate fluxes for 
the 140 observed H 2 emission lines and to reconstruct the Lya emission profile 
incident upon the warm H 2 . The excitation of H 2 can be determined from relative 
line strengths by measuring self-absorption in lines with low-energy lower levels, 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract NAS5-26555. This work is also based on data obtained for the FUSE Science Team by the 
NASA-CNES-CSA FUSE mission, operated by Johns Hopkins University. 



-2- 



or by reconstructing the Lya profile incident upon the warm H 2 using the total 
flux from a single upper level and the opacity in the pumping transition. Based 
on those diagnostics, we estimate that the warm disk surface has a column den- 
sity of log iV(H 2 ) = 18. 5^0 8) a temperature T = 2500^^ K, and a filling factor 
of H2, as seen by the source of Lya emission, of 0.25 ± 0.08 (all 2a error bars). 
TW Hya produces approximately 10~ 3 L Q in the FUV, about 85% of which is in 
the Lya emission line. From the H I absorption observed in the Lya emission, 
we infer that dust extinction in our line of sight to TW Hya is negligible. 

Subject headings: accretion, accretion disks — circumstellar matter — line: 
identification — stars: individual (TW Hya) - - stars: pre-main sequence - 
ultraviolet: stars 



1. INTRODUCTION 

Classical T Tauri Stars (CTTSs) are roughly solar-mass pre-main sequence (PMS) stars 
that are accreting gas from their circumstellar disks. The extent and mass of dust in disks 
have been determined from IR spectral energy distributions (SEDs) and imaging (see review 
by Zuckerman 2001). Cold gas at large radii is traced by observations of molecules such 
as CO and HCN (e.g., Aikawa et al. 2002). Warmer neutral and ionized gases are used as 
diagnostics of accretion (e.g., Gomez de Castro & Lamzin 1999; Johns-Krull & Gafford 2002). 
However, observing the gas in the disk close to the star, which is important for the formation 
and evolution of planets, has been difficult. Gas at these radii can induce planet migration 
(Goldreich & Tremaine 1979), dampen the eccentricity of terrestrial planets (Agnor & Ward 
2002), and is neccessary for accretion onto giant planets (Lissauer 1993). 

In search of gas at planet-forming radii, Calvet et al. (1991) investigated the formation 
of CO overtone lines at 2.3 fim in disks around CTTSs. These models predict that CTTSs 
as a group should show both CO emission and absorption, while Herbig Ae/Be stars should 
show CO in emission. Using high spectral resolution observations, Carr et al. (1993) showed 
that the CO overtone emission from the probable Herbig Ae/Be star WL 16 arises in the 
inner disk, confirming the prediction of Calvet et al. (1991). Najita et al. (1996) find the 
same result for the Herbig Ae/Be star 1548C27. However, broader searches at high spectral 
resolution for CO overtone emission or absorption from CTTSs generally find no evidence 
for contributions to these lines from the disks around these stars (Casali & Eiroa 1996; 
Johns-Krull & Valenti 2001a). On the other hand, Najita, Carr, & Mathieu (2003) find that 
CO fundamental emission near 4.6 - 4.9 /iin is common among CTTSs. They suggest that 
overtone emission has not been detected at the observed CO temperatures (1100 - 1300 K), 
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mainly because the overtone bands have much lower transition rates. However, CO is not 
expected to be the dominant gas constituent in the disks around CTTSs, and it is the major 
component of the gas, H 2 , that we seek to probe here. 

Previous observations of H 2 have not provided a clear picture of the gas in the disk. Thi 
et al. (2001) detected cold H 2 emission in the pure rotational S(0) and S(l) lines from ISO 
observations of a large sample of CTTSs, Herbig Ae/Be stars, and debris-disk stars. However, 
in ground-based observations of the S(l) and S(2) lines, using much smaller apertures than 
ISO, Richter et al. (2001) and Sheret et al. (2003) did not detect any H 2 emission from some 
of the same sources, indicating that the emission is probably extended. On the other hand, 
the emission in the H 2 1-0 S(l) transition detected by Bary, Weintraub, & Kastner (2003) 
in Phoenix spectra of 4 T Tauri stars is narrow and centered at the radial velocity of the 
stars, presumably formed in low-velocity circumstellar gas. 

While H 2 lines in the IR are typically weak and are observed against a strong dust 
continuum, fluorescent H 2 lines dominate the far-ultraviolet (FUV) spectrum of CTTSs at 
wavelengths longward of Lya. Ardila et al. (2002a) found that narrow H 2 lines can be 
blueshifted by up to 20 km s" 1 in HST /GRRS spectra of CTTSs, observed with a 2" x 2" 
aperture. High-resolution echelle and long-slit spectroscopy with #ST/STIS is providing the 
critical data needed to identify the source of the H 2 emission. Observations of T Tau reveal 
two components of H 2 emission. Saucedo et al. (2003) and Walter et al. (2003) detect off- 
source UV H 2 emission, that is pumped near line-center of Lya, and most likely produced 
where stellar outflows shock the surrounding ambient molecular material. The same set 
of fluorescent H 2 lines are detected in low-excitation HH objects such as HH43 and HH47 
(Schwartz 1983; Curiel et al. 1995), suggestive of a similar excitation mechanism. Walter 
et al. (2003) find that the on-source H 2 emission is photoexcited by a much broader Lya 
emission line, likely produced by the accretion shock at the surface of T Tau. 

In Herczeg et al. (2002, hereafter Paper I), we presented observations of H 2 fluores- 
cence in the UV spectrum of TW Hya obtained with HST/STIS and FUSE. Several lines 
of evidence suggested a disk origin for the H 2 emission: (i) the H 2 emission is not spatially 
extended beyond a point source in the cross-dispersion direction, (ii) the H 2 lines are emit- 
ted interior to TW Hya's wind because emission in one H 2 line is supressed by C II wind 
absorption, (iii) the H 2 line centroids have the same radial velocity as the photospheric lines 
of TW Hya, (iv) no H 2 absorption is detected against the Lya emission, and (v) the TW 
Hya Association is isolated from large reservoirs of interstellar molecular material, making 
a circumstellar origin for the H 2 emission unlikely. Because the H 2 emission is not extended 
beyond a point source, the spatial resolution of HST/STIS restricts the emitting region to 
be within 2 AU of the star. 
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Analysis of on-source IR H 2 emission lines to date has been limited because typically only 
one or two lines are observed. The IR lines could be excited by shocks, UV fluorescence, 
X-rays, or thermal heating. In contrast, over 140 H 2 lines from 19 different upper levels 
are detected in the FUV spectrum of TW Hya. The excitation of H 2 can be determined 
from relative line strengths in this rich spectrum either by measuring self-absorption in low 
excitation members of individual fluorescent progressions or by using an assumed Lya profile 
to reconstruct initial level populations in various pumping transitions. 

In this paper, we calculate the temperature of warm gas at planet-forming distances by 
modelling warm H2 in a disk around TW Hya. We use the H2 emission to reconstruct the 
Lya profile incident on the fluoresced H 2 . The observed Lya profile differs from the incident 
profile because of interstellar absorption along our line of sight to TW Hya. We discuss the 
strength of Lya emission and its affect on the circumstellar disk. 

2. PROPERTIES OF TW HYA 

TW Hya is the namesake of the 10 Myr old TW Hya Association (TWA), which most 
likely originated in the Sco-Cen OB association (Mamajek, Lawson, & Feigelson 2000). TW 
Hya is one of the oldest known stars still in the CTTS phase, and its mass accretion rate of 
about 2 — 3 x 10~ 9 M Q yr -1 is lower than the mass accretion rate of most of the 1-2 Myr 
old CTTSs that are found in Taurus (Kenyon & Hartmann 1995). From PMS evolutionary 
tracks, Webb et al. (1999) estimate that TW Hya has a mass of 0.7 M & and a radius of 1.0 
R & . Gaseous material is often associated with young stars, but is not found in the nearby 
TWA, so the extinction is low. The TWA is located at a distance of approximately 56 pc 
(Webb et al. 1999), compared with 140 pc for Taurus. Consequently, TW Hya is by far the 
brightest known CTTS in the UV and in X-rays. Most of the UV radiation from TW Hya 
is presumably produced at the accretion shock. The disk is viewed face-on (e.g., Zuckerman 
et al. 1995; Alencar & Batalha 2002), which prevents signficant Keplerian broadening in the 
H 2 line profiles. The disk mass is roughly (1.5 - 3) x 10~ 2 M© (Wilner et al. 2000; Trilling 
et al. 2001) and extends more than 225 AU from the star (e.g. Krist et al. 2000). Based on 
estimates of the stellar mass, radius, and rotation period of TW Hya, Johns- Krull & Valenti 
(2001b) calculate that the disk at 6 i?* corotates with the stellar surface. In the Shu et al. 
(1994) model, the disk is truncated at corotation. Weinberger et al. (2002) estimate from 
modelling the IR SED of TW Hya that the dust in the disk is truncated at about 11 R*. 
In their models of the IR SED, Calvet et al. (2002) found an underabundance micron-sized 
dust grains within 4 AU of TW Hya. Figure 1 shows a schematic model of TW Hya and its 
disk. 
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Since the mass accretion rate of CTTSs can vary, we calculate the mass accretion rate 
at the epoch when we observed TW Hya. We then proceed to estimate the extinction using 
H I absorption against Lyct emission and the SED of the NUV continuum that is produced 
in an accretion shock. The observations discussed in this paper were described in Paper I 
and are listed in Table 1. 



2.1. Mass Accretion Rate 

The strong excess blue and UV emission from CTTSs can be modelled as emission from 
shocks at accretion footpoints on the star, which leads to estimates of the mass accretion 
rate. We estimate the accretion rate of TW Hya by following the method of Valenti, Basri, & 
Johns (1993), who modelled observed blue spectra with the sum of a photospheric template 
and emission from a slab of hot isothermal hydrogen to crudely simulate emission from the 
accretion shock. 

We use the HST/STIS spectrum of the K7 weak-lined TTS V819 Tau (see Table 1), 
spanning 1100-6000 A, as a non-accreting template for the photospheric emission from TW 
Hya. For V819 Tau we adopt a visual extinction of Ay = 1.52, which is the average of values 
determined by Kenyon & Hartmann (1995) and White & Ghez (2001). We then use the 
interstellar extinction law of Cardelli et al. (1989) to remove the effects of extinction from 
the observed spectrum of V819 Tau. The NUV and FUV emission from V819 Tau is much 
fainter than that of TW Hya, and is not significant in this analysis. 

The spectrum of the model slab is determined by the slab temperature T, density n, 
thickness I, extinction Ay towards the slab, and a slab surface area parameter 5, as described 
in Valenti et al. (1993). The slab surface area parameter 5 corresponds to the fraction of 
the stellar surface area covered by accretion-related emission, and is calculated from the 
scale factor applied to the slab spectrum. The scaling factor between the non-accreting 
template star and TW Hya is another free parameter. We find the best-fit model parameters 
by fitting the synthetic slab spectrum plus the scaled template spectrum to the NUV and 
optical spectrum of TW Hya, minimizing x 2 using the amoeba function in IDL, which is 
based on the downhill simplex minimization method. The best-fit parameters for the slab 
are A v = 0.0, T = 9900 K, n = 1.7 x 10 14 cm" 3 , / = 125 km, 5 = 0.008, and a template 
scaling factor of 1.18. 

Figure 2 shows our fit to the optical spectrum of TW Hya, after subtracting the scaled 
template spectrum of V819 Tau. The accretion continuum calculated in this model is con- 
sistent with the observed NUV continuum. However, the observed continuum flux rises 
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shortward of 1700 A, which could result from (i) an H 2 dissociation continuum produced by 
electron collisions (Liu & Dalgarno 1996), (ii) H 2 fluorescence pumped by the FUV contin- 
uum, or (iii) by a hot accretion or activity component. The SED of the FUV continuum is 
similar to the continuum observed in HH2, which is probably produced by H 2 dissociation 
(Raymond et al. 1997). However, the FUV continuum of TW Hya appears smooth and does 
not show strong H 2 emission lines at 1054 A and 1101 A, which are detected towards HH2. 

Assuming that half of the potential energy of the accreting gas is converted into hydrogen 
continuum emission 2 , then the mass accretion rate of TW Hya can be estimated by 

• 1.25i2*L B i ab 

where L s i ab = 3.5 x 10~ 10 erg cm" 2 s" 1 is the slab flux integrated over all wavelengths. The 
leading factor of 1.25, derived by Gullbring et al. (1998) for a magnetospheric accretion 
geometry, replaces the factor of 2 from a spherical accretion geometry. We adopt 0.7 M Q 
and 1.0 R Q for the mass and radius of TW Hya (Webb et al. 1999). Uncertainties in the 
stellar mass, radius, and distance are all of order 20%. We calculate an accretion rate onto 
the star of M ~ 2 x 10 -9 M & yr" 1 , which is consistent with the accretion rate of 2 x 10~ 9 M Q 
yr _1 estimated by Alencar & Batalha (2002) using a similar method, and is larger than the 
accretion rate of 5 x 10~ 10 M Q yr -1 estimated by Muzerolle et al. (2000) using an analysis 
of the Ha line. 



2.2. Extinction to TW Hya 

In this section, we argue that the extinction to TW Hya is negligible by measuring the 
hydrogen column density directly from the shape of the observed Lya line, and then using 
an ISM extinction law to infer the dust extinction. This argument does not necessarily imply 
that the extinction towards the warm H 2 is negligible, because some dust could be mixed 
with the warm H 2 , leading to extinction of H 2 fluxes without any extinction of the stellar or 
accretion emission. We will discuss this possibility in §5. 

Figure 3 shows the observed Lya profile, which is likely produced in an accretion column, 
contains a dark, broad absorption feature that extends from +200 km s _1 to -500 km s _1 . 
The absorption is produced by H I in the wind and in interstellar and circumstellar material. 
Because circumstellar and interstellar absorption are indistinguishable in our data, we model 
these two features with a single absorption feature. All UV spectra of CTTSs show similar 



2 See Lynden-Bell & Pringle (1974) and Hartigan et al. (1991) for details. 
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wind absorption in Mg II (Ardila et al. 2002b), so the star must be occulted by the outflow, 
regardless of inclination angle. 

In the following analysis, we fit two Gaussian emission profiles and either one or two 
Voigt absorption profiles to the observed Lja profile. The intrinsic emission profile is as- 
sumed to be a double- Gaussian profile, and the parameters of the Gaussians are allowed to 
vary to best fit the observed profile when combined with the absorption. Based on the widely 
separated wind and interstellar absorption features observed in the O I 1302 A, C II 1334 A, 
and Mg II 2795 A lines, the interstellar and wind absorption components towards TW Hya 
must be centered at v — and —130 km s" 1 , respectively. For the interstellar component 
we assume a Doppler parameter of b = 11 km s _1 , which corresponds to a temperature of 
7000 K that is typical of the local (warm) ISM (Redfield & Linsky 2000). An interstellar 
column density is then determined by fitting the absorption longward of line center, where 
any wind absorption should be negligible. The wind absorption extends to a velocity of at 
least -500 km s" 1 , and likely has a large Doppler parameter due to velocity gradients along 
our line of sight. Given these assumptions, the interstellar neutral hydrogen column density 
in our line of sight is log A(H I) = 19.05 ± 0.15 (Fig. 3). The upper limit of the column 
density in the wind is A(H I) < 19.4 for Doppler parameters b > 25 km s _1 

In a more conservative test to find the upper limit to the hydrogen column density, we 
fit the Lya absorption with a single component, effectively combining any interstellar and 
wind absorption. We vary the central wavelength of the absorption and set a low Doppler 
parameter parameter, so that the damping wings generate the width of the absorption pro- 
file. Figure 4 shows that the maximum possible column density consistent with the data is 
logA(H I) = 19.75. 

The conservative upper limit of the hydrogen column density towards TW Hya is there- 
fore log AT (H) < 19.75. In the only previous measurement of N(H I) to TW Hya, Kastner 
et al. (1999) used X-ray spectral energy distributions from ROSAT and ASCA to estimate 
higher column densities of logA"(H I) = 20.5 ± 0.4 and 21.4 ± 0.1, respectively. How- 
ever, the X-ray analyses are highly reliant on the accuracy of plasma codes. Kastner et al. 
(1999) estimate the hydrogen column density by assuming solar abundances and using a 
two-temperature emission measure rather than a differential emission measure distribution. 
However, in their analysis of the Chandra X-ray spectrum of TW Hya, Kastner et al. (2002) 
found that the abundances show an inverse-FIP effect. We further note that the sensitivity 
and flux calibration of ROSAT and ASCA are relatively poor at soft energies, where ab- 
sorption from H I ionization is measured. Figure 4 shows that the absorption calculated for 
the Lya line for log A(H I) = 20.5 is signficantly broader than the observed absorption pro- 
file. Our comparison of the calculated transmission function with the observed Lya profile 
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provides a much simpler and more direct test of hydrogen column density than the X-ray 
analyses. The discrepancy between the UV and X-ray measurements of N(H) could in princi- 
ple be due to a variable circumstellar absorber, which Kastner et al. (1999) invoked to explain 
the order of magnitude discrepancy in calculations of N(H) from ASCA and ROSAT. How- 
ever, the many ultraviolet observations of TW Hya by IUE, HST '/STIS, and FUSE, which 
are widely spaced in time, show no evidence for large fluctuations in extinction. 

We estimate an H 2 column density in our line of sight to TW Hya with FUSE, as the 
bandpass contains many H 2 transitions from low rotational levels of the ground electronic 
state, which are often seen in absorption due to the ISM at temperatures typically ~ 100 K 
(Rachford et al. 2002). H 2 absorption lines from v" = 0, J" = 0, 1, 2, 3, 5 can in principle be 
observed against O VI emission (e.g. Roberge et al. 2001). The O VI profiles are noisy (see 
Fig. 5), but we estimate an upper limit of log V(H 2 ) < 18 by calculating the transmission 
percentage for a range of temperatures and column densities, meaning the molecular fraction 
of H towards TW Hya is very low. 

When we convert total hydrogen column density, V(H I)+2V(H 2 ) = V(H), to a red- 
dening using the interstellar relationship of Bohlin, Savage, & Drake (1978), 

MHI + H 2 ) 9i 9 

E(B - V) = 5 ' 8 X 10 Cm (2) 

we find that the reddening toward TW Hya is E(B — V) = Ay / R v < 0.01, where Ry is 
a dimensionless parameter that ranges between 2.5-6, depending on the size distribution of 
dust grains. Significant deviations from the assumed dust-to-gas ratio can occur in star- 
forming regions such as p Oph, in which the dust grains are likely larger than they are in the 
ISM. However, the low molecular fraction towards TW Hya is consistent with a low extinction 
(Savage et al. 1977; Rachford et al. 2002), and is not consistent with material found towards 
star-forming regions. We conclude that this relationship provides a reasonable estimate of 
the extinction along our sightline to TW Hya. Even if the dust-to-gas ratio along the line 
of sight to TW Hya were for some reason enhanced by an order of magnitude, then dust 
extinction would reduce the brightness of TW Hya by a factor of 2 and have only a 20% 
effect on the relative fluxes between 1200 A and 1650 A. 

Inspection of the NUV continuum (Fig. 6) confirms the small extinction to TW Hya. 
If we assume that the NUV continuum is produced in an isothermal slab, as modelled in 
§2.1, then it can be well fitted with a hydrogen continuum only for extinctions Ay < 0.2. 
As noted in §2.1, an accretion continuum best fits the NUV continuum for Ay = 0.0. 
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3. MODELLING H 2 LINE FLUXES 

In Paper I, we measured a total line emission flux of 1.94 x 1CT 12 erg s -1 cm -2 by 
summing 146 H2 emission lines in the 19 progressions (the set of R and P transitions from 
a common upper level) listed in Table 2. Most of these lines are Lyman-band transitions 
fluoresced by the strong, broad Lja emission line. 

Wood, Karovska, & Raymond (2002) discovered in their analysis of HST/STIS obser- 
vations of Mira B that H 2 fluxes weaken towards short wavelengths because of increasing 
H 2 line opacities. This occurs because shorter wavelength lines are transitions to low en- 
ergy levels, which are more populated than higher levels. Figure 7 demonstrates that the 
same effect occurs in the H 2 fluorescence sequences detected towards TW Hya. In most 
sequences, fluxes in many H 2 emission lines at short wavelengths are weaker than predicted 
from theoretical branching ratios because of significant absorption in the same transitions. 
Transitions from levels with low excitation energies, such as v'—0, J'=l, have significant 
optical depth below 1300 A, but this does not happen for transitions from levels with large 
excitation energies such as v'—A, J' =18 (see Fig. 7). The v'=4, J' =18 level does not have 
any downward transitions with significant opacity in the STIS bandpass. 

Following the Monte Carlo method of Wood et al. (2002), we compute H 2 opacities by 
modelling H 2 fluorescence in a ID isothermal plane-parallel atmosphere of H 2 irradiated by 
Lja photons. We assume complete frequency and angular redistribution when photons scat- 
ter in the slab. Forty-one wavelength grid points across each H 2 line profile track differences 
in escape probability across the line profile because because a photon that scatters into the 
wings of a line escapes from a slab more easily than a photon at line center. A Lya photon 
enters the slab at an angle 9 with respect to the slab surface normal and photoexcites an H 2 
molecule. The H 2 molecule quickly de-excites to the ground electronic state through one of 
the many available routes, or dissociates with a probability Pois [listed in Table 2, calculated 
by Abgrall, Roueff, & Drira (2000)]. Depending on the optical depth of H 2 in the lower level 
of the downward transition, the resulting photon may be temporarily absorbed by another 
H 2 molecule or may emerge from either side of the slab at an angle <fi with respect to the 
normal. In our simple model, H 2 line photons created fluorescently in the slab either escape, 
are converted to another H 2 line in the same progression, or dissociate an H 2 molecule. No 
other photon processes are considered. We tabulate photons that exit the slab on the Lja 
entry side (reflected photons) and opposite side (transmitted photons) using 40 bins of equal 
solid angle on each side of the slab. We refer the reader to Wood et al. (2002) for further 
details of this Monte Carlo code. 

We use the term "model" to refer to the simulated plane-parallel atmosphere described 
in the preceding paragraph, and the term "geometry" to refer to the morphology of the H 2 
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around the Lja emission source. In the context of the model, the H 2 line fluxes depend on 
the temperature and column density of the slab, the angle of incidence of the Lja photon 
into the disk, the outward angle of the emerging photon, the side of the slab from which 
the photon escapes, and the optical depth of the transition. We assume that extinction is 
negligible in the slab and in any interstellar material. Figure 8 shows two possible geometries 
that we use to approximate a thick disk and a thin disk. The term "thick disk" applies to 
a star-disk geometry in which most of the H 2 illuminated by Lja is at the inner edge of 
the disk. The term "thin disk" applies to a geometry in which most of the H 2 illuminated 
by Lja is at the disk surface. In this section and §4, we present detailed results for the 
thick disk geometry that reproduces the data well. In this adopted geometry, Lja photons 
enter the slab normally (6 = 0°), and the H 2 photons that we observe emerge at an angle 
= 76° — 79° to the normal. This choice is intended to present one approximation of a 
plausible disk geometry for which we can present detailed results, and is not meant to exclude 
other geometries. In §5, we present general results for many geometries to demonstrate that 
our primary result, the temperature of the H 2 gas, is relatively insensitive to geometries, 
including those that are significantly different from the adopted geometry. 

We calculate the relative flux of each transition within a single progression, given an H 2 
excitation temperature and column density of the slab, and a separate flux normalization 
factor for each progression. Models are calculated for temperatures in the range T = 1000- 
5500 K, sampled at 100 K intervals, and column densities (in cgs units) of log 7V(H 2 )=15-23, 
sampled at 0.1 dex intervals. iV(H 2 ) and T are then determined by fitting simultaneously 
the relative line fluxes for 19 observed progressions. Five lines are not used because they 
are blended, and one line is not used because it is anomalously weak due to C II wind 
absorption (see Paper I). The FUSE lines, observed six weeks after the STIS observations, 
are not used in building these models due to the potential problems with flux variability 
and cross-calibration between STIS and FUSE. The fits have 21 variables: 19 normalization 
factors, iV(H 2 ), and T. In typical models, 46 of the 140 H 2 line fluxes are sensitive to H 2 line 
opacity, and hence iV(H 2 ) and T. The other 94 lines, which have only a weak dependence 
on iV(H 2 ) and T, determine the 19 normalization factors. We calculate Xm f° r on ly the 46 
lines which depend on iV(H 2 ) and T. Confidence contours of 1 — 5a, calculated from values 
of Axl m for 44 degrees of freedom, are shown in Figure 9, with a minimum xtm = 2-22 at 
logiV(H 2 ) = 17.9 and T = 2600 K. 

Based on this analysis, and additional constraints provided by the Lya profile analysis 
described in §4, we ultimately decide a thick disk model with T = 2500 K and logiV(H 2 ) = 
18.5 is our best model, which we adopt for the analysis (see §6.1). Figure 10 compares 
the observed spectrum compared with a synthetic spectrum computed from this model. The 
adopted model correctly predicts the fluxes of all strong H 2 lines. Only one observed line flux 
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differs significantly from the model flux. The model predicts a flux of 46 x 10 15 erg cm 2 s 1 
in the 0-4 R(l) transition at 1333.797 A. The measured flux in this line, (7.9 ±0.7) x 10~ 15 
erg cm -2 s -1 , is significantly reduced by C II wind absorption, as discussed in Paper I. 
The 0-3 P(3) transition at 1283.16 A, with a measured flux of 28.0 x 10~ 15 erg cm~ 2 s _1 , 
is moderately weaker than the model flux of 34.8 x 0" 15 erg cm~ 2 s _1 , possibly because of 
H 2 absorption in the 2-2 R(15) transition, centered 8 km s _1 from the 0-3 P(3) transition. 
Table 3 lists predicted fluxes for Lyman-band H 2 transitions in the FUSE bandpass. Future 
measurements of flux in these transitions would provide a good test for our model. Lines at 
A < 1200 provide particularly good constraints since they are transitions to very low energy 
levels, and as a result have large opacities. 

Neither this model nor the alternative models presented in §5 rigorously simulate a 
disk geometry. Because the structure of disks is uncertain near the truncation radii, we 
do not have an accurate picture of the geometry of the H 2 emission region. Without this 
a priori knowledge, the thick disk geometry described here is not clearly better or worse 
than any other oversimplified geometry that we might assume. In §5, we will show that 
the temperature of the molecular region is relatively insensitive to geometry, and we will 
rule out certain geometries that cannot sufficiently explain the observed H 2 spectrum. Since 
our models make as few assumptions as possible, the primary result of this paper, the 
temperature of the H2 gas, may apply to a wide range of other morphologies. 



4. RECONSTRUCTED Lya PROFILE 

In this section, we reconstruct the Lja profile seen by the warm H 2 gas at the 18 
pumping wavelengths between 1212-1220 A. Fluorescent H 2 emission occurs because certain 
Lyman-band transitions, with lower levels 1 — 2 eV above the ground state, have wavelength 
coincidences with the broad Lja profile. The emission in each progression depends on the 
Lja flux at the pumping wavelength of the fluorescent transition, the upward transition 
probability, the opacity in the pumping transition, and the filling factor of H 2 as seen from 
the source of the Lja emission. We reconstruct the incident Lja profile following a pro- 
cedure similar to that used by Wood et al. (2002) to analyze fluorescent H 2 emission in 
the HST /STIS spectrum of Mira B. A comparison of the reconstructed and observed Lja 
profiles provide another constraint on the values of T and iV(H 2 ), that is independent of the 
constraints provided by the analysis presented in §3. The analysis in §3 relied upon using 
the H 2 flux ratios within each fluorescence progression to constrain T and iV(H 2 ), whereas 
the Lja profile reconstruction analysis presented in this section relies only on the total flux 
from each upper level. 
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We calculate the optical depth in each H 2 pumping transition by assuming that the 
absorption line profiles in the model slab are Voigt functions with a Doppler b parameter, 
which depends on an assumed temperature, and a damping parameter a from Abgrall et al. 
(1993). The total Lya flux, F abs , absorbed by a given lower level is then 

F abs = V x F inc x J 1 - e-^PW^dA (3) 

where F inc is the incident Lya flux per unit wavelength and f x 1 — exp" 

-r x [T,N(H2)] dX j s the 

equivalent width of the absorption profile within the slab (see EQW^ in Table 2). Not 
every absorbed Lya photon produces an observed H 2 photon. For the most highly excited 
states, the dissociation probability (P dis in Table 2) can be as large as 42%. Some fraction 
(Runseen in Table 2) of the emitted photons will be in spectral lines that we did not observe 
due to line blends, noisy data, or insufficient wavelength coverage. In our plane-parallel slab 
model, an emerging photon travelling normal to the slab sees a smaller optical depth than 
a photon travelling nearly parallel to the slab. As a result, the emitted flux in these models 
is largest for angles normal to the slab. 

For comparison with the observed Lja profile, we compute reconstructed Ly-alpha 
profiles for the same grid of T and iV(H 2 ) values used in §3. Model fluxes at 8 pumping 
transitions, all on the red wing of the Lya line, are fitted to the observed profile. The observed 
Lja flux at other pumping wavelengths is corrupted by interstellar or wind absorption. Two 
progressions, 3-3 R(2) at 1217.031 A and 3-3 P(l) at 1217.038 A, are not included in the 
fit because their absorption profiles overlap, complicating the analysis. Photoexcitation of 
H 2 by C III emission via the transition 3-2 P(4) at 1174.923 A may also add flux to the 
progression pumped by 3-3 R(2). We judge the suitability of a model using a modified x 2 
analysis to measure the scatter in the red wing of the reconstructed profile. Errors in the 
Lja line wing tend to be multiplicative, rather than additive, so we define x\ va 




based on standard error propagation techniques from Bevington (1969). Using this formula, 
an overestimate and an underestimate of the reconstructed Lya profile by a factor of 10 will 
result in the same contribution to xl ya - 

Fi gure 11 shows confidence contours of 1 — 5cr calculated from ^XuLya for fits to the 
observed Lja profile, with 5 degrees of freedom. In the thick disk geometry, a minimum 
Xu,Lya = 2-31 is obtained when comparing the reconstructed and observed Lja profiles for 
log7V(H 2 ) = 18.5, T = 2500 K, and rj = 0.25. The top panel of Fig. 12 shows that these 
parameters yield an excellent fit to the observed profile. These parameters also reproduce the 
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individual line fluxes well (see §3 and Fig. 9). Since the confidence contours calculated from 
the Lya profile analysis place a more stringent constraint on the temperature and column 
density than do the contours from fitting individual line fluxes (see Fig. 11), we choose on 
the values of T and 7V(H 2 ) quoted above as our best estimates of these parameters. The 
model dependent quantities listed in Table 2 are all calculated using these parameters. 

More than 500 other possible upward transitions in the Lyman-band have wavelengths 
between 1210-1220 A. Fluorescent emission in lines from these weaker progressions is below 
our detection limit. Table 4 identifies undetected progressions with the largest predicted 
emission line fluxes. The table includes our observational flux limit F max on the strongest 
line in the progression at wavelength \ ma x- For any H 2 line that is not detected and does 
not appear in Table 4, we place an upper flux limit of 5 x 10~ 15 erg cm -2 s -1 . We then use 
these upper limits to calculate an upper limit of Lya at the pumping wavelength to test the 
self-consistency of our model. 

Figure 12 shows three examples of reconstructed Lya profiles (circles), with upper limits 
to the Lya flux at certain wavelengths calculated from upper limits of lines in the undetected 
progressions (arrows), as described above. The top panel shows the Lya profile constructed 
using the parameters for our adopted model [log jV(H 2 )=18.5, T = 2500 K, and rj = 0.25]. 
This profile is smooth, and its wings correspond well with the wings of the observed Lya 
profile. The upper limits on the flux of Lya calculated from undetected progressions are 
all consistent with the Lya profile reconstructed from observed progressions. The middle 
panel shows a Lya profile calculated from a model with log 7V(H 2 )=19.0, T = 2000 K, and 
7] = 0.4, which is inconsistent with the data because of a large amount of scatter in the red 
wing of the reconstructed Lya profile. The bottom plot shows a Lya profile calculated from 
a model with log 7V(H 2 )=19.0, T = 3000 K, and rj = 0.2, which is also inconsistent with the 
data because at 1214.421 A the predicted Lya emission required to pump the v' — 7,J' — A 
level is significantly below the observationally constrained flux at 1214.465 A. Figure 13 
is a reproduction of Figure 11, with dotted and diagonal lines added to show the region 
of parameter space that are excluded because upper limits on H 2 emission would imply a 
Lya pumping flux significantly lower than the nearby Lya pumping flux for a progression 
of observed H 2 lines. Parameter space constraints are shown separately for upper limits on 
progressions pumped via 7-4 P(5) at 1214.421 A alone (dots) and via a set of others upward 
transitions (diagonal lines). The parameter space with r\ > 1 (dashed lines in Fig. 13) are 
unphysical models, and are also ruled out. 

The three reconstructed Lya profiles shown in Figure 12 are roughly similar to each 
other and to most other Lya profiles that we have reconstructed for other values of T and 
jV(H 2 ). Each profile shows a clear depression near Lya line center that is not nearly as 
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wide as the observed H I absorption. Figure 14 demonstrates that the absorption feature in 
the reconstructed Lya profile can be characterized by a column density logiV(H I) m 18.7 
centered 90 km s^ 1 shortward of line center. This H I between the source of Lya emissio and 
the warm H 2 could be in the base of a wind or at the disk surface. Alternatively, the shape 
of the Lya profile could be an intrinsic self- reversal, in which case the intervening H I need 
not exist. The blue side of the reconstructed Lya profile is enhanced relative to the observed 
Lya profile, because wind absorption in our line of sight to the star attenuates the intrinsic 
emission. This result supports our conclusion from Paper I that, in our line of sight, the H 2 
must be photoexcited interior to most of the wind that we detect. 



5. ALTERNATIVE GEOMETRIES 

In §3 and §4, we described in detail results for a dust-free thick disk geometry, which 
may not represent the true morphology of warm gas around the star. In this section, we 
explore the same parameter space for a thin disk geometry, shown schematically in Figure 
8, and a cloud geometry, which approximates H 2 emission from a cloud in our line of sight 
to the star. 

The computed H 2 fluxes depend on the angle at which Lya photons enter the disk, the 
exit angle of H 2 photons, and the extinction either within the disk or in the ISM. Figures 
15 and 16 and Table 5 present the best-fit parameters from fitting the H 2 line opacities for 
each of three geometries as a function of the exit angle of H 2 and extinction within the disk. 
In general, while changes in the geometry of the system can mildly change some results, the 
calculated temperature of the warm gas depends primarily upon the molecular physics and 
not on the geometry. The calculated column density of H 2 and its physical interpretation, 
however, do depend on the assumed geometry. In what follows, the column density represents 
the amount of H 2 in a line of sight normal to the disk surface, regardless of the path length 
the Lya photons traverse through the disk. 



5.1. Thin Disk Geometry 

A thin disk (see Fig. 8) is characterized here as a disk for which the filling factor of H 2 
around the Lya source is dominated by the surface layers of the disk, rather than by the 
inner edge of the disk. We approximate this geometry with a model in which the incident 
Lya photons enter a plane-parallel slab at a large angle with respect to the disk normal 
and are reprocessed into H 2 photons that emerge normal to the slab. Because Lya photons 
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enter the slab at a large angle, the radiation is reprocessed into H 2 emission at a shallow 
depth, and the emerging H 2 photons see a smaller optical depth in this geometry than in 
the thick disk geometry. For the thin disk geometry, with incident photons entering the disk 
at — 80°, the best-fit parameters are T = 2500 K, iV(H 2 )=20.0, but with an unacceptable 
xIh2 = 7-9, for fits to H 2 line fluxes. For fits of the reconstructed to the observed Lja line, 
the confidence contours move to smaller column densities by a factor of sin 9. As the angle 
of incidence asymptotically approaches 90°, the Lja photons are converted into H 2 photons 
close to the disk surface. Since the optical depth traversed by the emitted H 2 photons is 
negligible, the relative line fluxes with each progression differ from predicted branching ratios 
only because of extinction, which we find to be unlikely (see §2.2 and §5.3). 



5.2. Cloud Geometry 

We next consider a geometry in which the H 2 is located between us and the source of 
Lja photons. In this geometry, the escaping photons must travel through a greater optical 
depth of H 2 than they would in the other geometries that we considered. As a result, the 
acceptable contours for fits to H 2 line fluxes move towards smaller column densities relative 
to the thick or thin disk geometries. Even in this very different geometry, the confidence 
contours for the Lja reconstruction fitting do not move much in parameter space, and they 
restrict the gas temperature to 2000-3000 K. The best-fit parameters in fitting H 2 fluxes for 
this geometry are T = 2600 K, iV(H 2 ) = 17.8, with xl,H2 = 2-8. The transmission geometry 
is probably not realistic for TW Hya, but it can be used to set limits on the column of warm 
H 2 between us and the source of Lya emission. The absence of any observed H 2 absorption 
in the Lya profile (Fig. 3) allows us to place a limit on the column density in this level of 
log N(v" = 2, J" = 1) < 14.4 and a total column density of warm H 2 of log iV(H 2 ) < 16.1. 
Thus, the detected H 2 fluorescence does not occur in our line of sight to TW Hya. 



5.3. Dust Extinction in the Disk and the ISM 

In §2, we presented evidence that the dust extinction along our line of sight to the Lja 
emission source is negligible. Extinction of H 2 emission could occur within the disk if dust 
is mixed with the warm gas, although Calvet et al. (2002) found that small dust grains 
are underabundant within 5 AU of the star. We assume an interstellar extinction law from 
Cardelli et al. (1989), for grains typical of the interstellar medium (Ry = 3.09). 

For each H 2 progression, the average depth d a b s of the initial Lja absorption depends 
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on the opacity of the lower level of H 2 (see Fig. 17 and Table 2). If a Lya photon is, on 
average, absorbed close to the disk surface, then the H 2 lines excited by that photon will 
not be significantly attenuated by any disk extinction. Thus, the extinction varies for the 
different progressions. In Figure 16, Ay represents the extinction through half of the slab 
(not the extinction through the entire slab or the extinction from the star to the observer), 
which is equal to the largest possible average extinction of H 2 emission when Lya photons are 
absorbed uniformly throughout the slab. Disk extinction attenuates both the incident Lya 
emission and the outgoing H 2 emission, but not the observed Lya emission, and thus increases 
the already large filling factor. In the context of our model, significant disk extinction can 
be ruled out because the computed filling factors would be much larger than unity. 

Interstellar extinction affects the fluxes smoothly across the wavelength region, and does 
not significantly impact the derived filling factor because it attenuates both the observed Lya 
emission and the H 2 emission. If we consider only extinction and ignore opacity in the H 2 
lines, then from fitting H 2 line fluxes we obtain a minimum xIh2 = 6-06 at Ay =1-3, which 
is significantly larger than the xl,H2 f° r the best-fitting models. 

5.4. Adoption of the Thick Disk Geometry 

Based on these results, and the results presented in Figs. 15-16 and Table 5, we adopt 
the thick disk geometry (Geometry 2 in Table 5) as the preferred geometry to describe the 
H 2 emission. The thick disk geometry is physically plausible and produces reasonably good 
fits with self-consistent results, although certain other geometries can produce lower values 
of x 2 m fits to H 2 fluxes or in the Lya fits. Some of those geometries (e.g., Geomteries 7 
and 8 from Table 5) represent morphologies that are difficult to reconcile with our previous 
conclusion that the H 2 emission occurs in a disk. In other geometries (e.g., Geometries 1 
and 11 from Table 5), contours from the H 2 flux fitting do not overlap with the contours 
from the Lya fitting as well as they do for the thick disk geometry, i.e. the joint probability 
is lower (see Fig. 11). 

We caution the reader that the thick disk geometry is best only in the context of our 
restricted set of models. Warm H 2 in the TW Hya star-disk system could in principle have 
a morphology not considered here. 
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6. DISCUSSION 

6.1. Synthesizing the H 2 and Lya modelling results 

Figures 11 and 13, and Table 5, compare results from our analysis of individual H 2 line 
fluxes with our results from the Lya profile reconstruction. Based on the overlap region of 
acceptable models for both the individual H 2 line fluxes and the reconstructed Lya profile, 
we conclude that a column density of logAf(H 2 ) = 18.5lgg is heated to a temperature of 
25001^ K, with a filling factor of rj = 0.25 ± 0.08 around TW Hya (2a error bars). Table 2 
lists the percentage P^s of H 2 that dissociates from the upper level, the H 2 dissociation rate 
in terms of mass (M^ s ), and the total flux for each progression (F mo d) for our best-fit model 
of logiV(H 2 ) = 18.5 and T = 2500 K. The notion of a single temperature is a simplistic 
assumption, but it is sufficient to reproduce our data. The temperature derived here is 
warmer than the ~ 1100 — 1300 K temperatures seen in the fundamental CO emission by 
Najita et al. (2003). 

The mass column corresponding to log iV(H 2 ) = 18.5 is 10~ 5 g cm -2 , which is about 10~ 7 
times smaller than the mass column predicted to be within 1 AU by D'Alessio et al. (1999). 
This suggests that H 2 fluorescence probably occurs in a very shallow surface layer of the disk, 
which could be a warm, extended disk atmosphere. The Lya fluorescence process cannot be 
used to detect cold H 2 , which most likely represents the bulk of the mass in the disk. 

In our models, we assume that H 2 gas is thermalized, so that relative level populations 
depend only on temperature. This is reasonable because the gas density at the disk surface 
should be between njj = 10 6 — 10 8 cm -3 (Gomez & D'Alessio 2000), which is far larger than 
the critical density of n#=10 4 cm -3 at which collisions typically dominate over quadrupole 
radiative de-excitation. Below the critical density, collisions cannot repopulate excited H 2 
levels as fast as they radiatively deexcite. Fluorescent progressions pumped out of higher 
energy states would be weaker than expected from our thermal models, generating scatter 
that we do not see in our reconstructed Lya profile. This conclusion confirms the assumption 
of thermally excited H 2 gas used by Black & van Dishoeck (1987) and Burton et al. (1990) 
to model H 2 fluorescence. 



6.1.1. Filling factor of H 2 around the Lya source 

Figure 13 shows that the filling factor of the H 2 around the Lya emitting region is 
77 = 0.25 ± 0.08. This factor represents the solid angle of the sky subtended by H 2 as 
seen from the source of Lya photons near TW Hya, although other geometrical effects may 
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contribute. Neutral hydrogen mixed with the warm H 2 could scatter Lya photons, increasing 
the effective path length of the Lya photons and thereby reducing the calculated filling factor. 
Any directional dependence of the Lya emission, and the percentage of Lya emission visible 
to us versus the percentage visible to the warm H 2 , will also affect the calculated filling 
factor. A total of 1-2% of the intrinsic Lya emission is reprocessed into H 2 emission (see 
Table 6 and §4), which places a lower limit of 0.01 on the filling factor. 

For any flared disk geometry, the filling factor can be converted into a disk height by 
comparing the surface areas of a cylinder and a sphere. For rj = 0.25, the implied geometrical 
disk height from the midplane to the surface at distance R from the center of the star is 
H=0.25(y^j) AU, compared with predicted disk heights from (D'Alessio 2001) of about 
0.15 AU at a distance of 1 AU from the star. 



6.1.2. Relationship between UV and IR H 2 Emission 

Emission in the 1-0 S(l) rovibrational transition at 2.1218 /zm has been detected from 
4 T Tauri stars (Bary et al. 2003). The narrow width and absence of any velocity shift with 
respect to the stars suggests a disk origin of this emission. The flux in this line from an 
optically thin emission region 

N(l,3)7ir 2 A ul he 

And 2 A ' [ ) 

where N(l,3) is the H 2 column density of the disk in the v" = 1, J" = 3 level, nr 2 is the 
projected surface area of the emission region, and A u i = 3.47 x 10~ 7 s _1 . If the 1-0 S(l) line 
flux of 1 x 10 -15 erg cm -2 s _1 measured by Weintraub et al. (2000) comes from a heated 
disk surface layer with N(H2)=18.5 and T=2500 K, then the 1-0 S(l) emission must extend 
out to at least 5 AU for a face-on disk. This corresponds to an angular extent of ±0'/l, 
which could have been resolved in our HST observations, but was not. Thus, it seems likely 
that at least some of the 1-0 S(l) emission comes from H 2 gas that is not fluoresced by Lya. 
Some of the emission in the 1-0 S(l) line may be produced deeper in the disk or farther 
from the star, as a result of X-ray ionization and subsequent collisional excitation of H 2 by 
non-thermal electrons (Weintraub et al. 2000). 



6.1.3. H 2 Dissociation Rate 

From the flux in observed transitions, we calculate a dissociation rate of 9.8 x 10~ n 
M Q yr _1 for H 2 in the disk. Forward modelling of undetected Lyman and Werner-band 
transitions excited by Lya photons increases the calculated dissociation rate to 1.6 x 10~ 10 M Q 



-19- 



yr -1 . Depending on the timescale for collisional de-excitation of H 2 , the actual dissociation 
rate could be much higher due to multiple pumping [pumping from lower rovibrational levels 
of H 2 that are populated by the initial pump and subsequent fluorescence - see Shull (1978) 
for details], which can populate highly excited rovibrational levels in the B electronic state 
that have dissociation probabilities as high as 50%. 

The H 2 dissociation rate we estimate is within a factor of 10 of TW Hya's mass accretion 
rate (see §2.1), meaning that the dissociation provided by Lya fluorescence may be an 
important process in the accretion of hydrogen onto the star. The neutral H released by 
dissociation will cause the incident Lya to scatter, increasing the probability of absorption 
in the pumping transitions and therefore decreasing the required filling factor. The neutral 
H could also be sufficiently optically thick to produce the self-reversal of the Lya profile that 
irradiates the disk (see Fig. 14). 

6.2. FUV radiation field 

The FUV radiation field incident on the accretion disk can significantly affect the chem- 
istry and temperature structure of the disk. In Table 6 we list the flux for spectral features in 
the FUV spectrum of TW Hya. The observed Lya line contributes 67% of the observed FUV 
flux of 6.5 x 10 -11 ergs cm -2 s _1 between 1170-1700 A. From the reconstructed Lya profile, 
we infer an intrinsic flux of 8 — 16 x 10 -11 ergs cm -2 s^ 1 , depending on the correction for 
the central absorption. With this correction, Lya emission contributes 80-90% of the FUV 
flux of TW Hya, of which 1-2% is reprocessed into H 2 emission. The continuum accounts for 
about 6 x 10 -12 ergs cm~ 2 s _1 , or about 4% of the total emission, including the estimated 
intrinsic Lya flux. The remaining flux occurs in other emission lines such as the C IV 1550 A 
doublet. The UV flux below 1170 A measured by FUSE adds at most 2% additional flux to 
the strength of the radiation field. The UV radiation between 1700-2000 A may add to the 
radiation field, primarily due to emission in lines such as Si II] A1817 and C III] A 1909, but 
will have little effect on the excitation of molecules and atoms in the disk. The estimated 
total emission after including Lya flux that is not observable is approximately 10~ 3 L Q , 
assuming isotropic emission and a distance of 56 pc. 

The properties of the disk, both in the PDR-like disk atmosphere and in the cooler outer 
regions, will be partially controlled by the Lya emission, as evidenced by the fluorescent H 2 
emission. Not only does Lya dominate the FUV flux of TW Hya, but it can also scatter 
off neutral H towards the interior of the disk, allowing it to penetrate the disk further than 
continuum UV radiation. Lya can dissociate molecules such as H 2 and H 2 0, and can ionize 
Si and C. Bergin et al. (2003) explain that the enhancement of CN relative to HCN occurs 
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in disks because Lya radiation can photodissociate HCN, whereas CN can be dissociated 
only by radiation with A < 1150 A. 

6.3. Heating mechanisms 

We have found evidence for a warm surface layer of H 2 on the disk around TW Hya. 
From the lack of spatial extent in the STIS echelle cross-dispersion direction, the heated 
surface is located within 2 AU of the star. Models of dust in disks of CTTSs typically 
predict T « 1000 K only within 0.1 AU of the star (D'Alessio et al. 2001; Chiang & Goldreich 
1997). The FUV radiation field at 1 AU from the star is approximately 10 7 Go, where the 
Habing field Go — 1-6 x 10~ 3 erg cm -2 s^ 1 (the local interstellar field is 1.7 Go), which is 
stronger than the FUV radiation field incident upon most PDRs (10 3 — 10 6 Go). Hollenbach, 
Yorke & Johnstone (2000) predicted from models of photo-dissociation regions that a warm 
atmosphere-like structure, heated by UV radiation, should surround the inner regions of a 
disk, with a temperature between 500 and 3000 K. 

Because of the absence of micron-sized dust grains within 5 AU of TW Hya (Calvet et 
al. 2002), a large percentage of the FUV radiation may be deposited into the gas. In their 
analysis of the H 2 fluorescence from Mira B, Wood et al. (2002) found that the temperature 
increases sharply with radius from the UV source, and that the Lya fluorescence process is 
a significant heating source for the gas. The energy deposited into the warm gas near TW 
Hya due to Lya photoexcitation and subsequent fluorescence is 1.4 x 10 29 erg s _1 , which 
corresponds to roughly 0.01 eV s™ 1 per molecule of H 2 . The strong X-ray luminosity of TW 
Hya, L x = 3.6 x 10~ 4 L Q between 0.45-6.0 keV (Kastner et al. 2002), could heat the disk 
surfaces by an additional 500 K (Igea & Glassgold 1999; Glassgold & Najita 2001; Fromang 
et al. 2002) by ionizing atoms to produce a large reservoir of energetic non-thermal electrons 
(e.g., Maloney et al. 1996; Yan & Dalgarno 1997). X-rays can also penetrate deeper into 
the gas than UV radiation and, consequently, may heat gas deeper in the disk than the UV 
radiation. 

Certain observed fluorescence progressions are pumped from lower rovibrational levels 
with large energies, such as v" = 3, J" = 25 (E" = 4.2 eV) and v" = 5, J" = 18 (E" = 3.8 
eV) . The measured flux in these progressions, comibed with the observed flux at the pumping 
wavelengths, requires that of order 1% of the H 2 resides in these levels. This population is 
too large to be thermally populated at any reasonable temperature for the H 2 gas. Therefore, 
non-thermal processes must be populating these levels. Although the fluorescence process can 
populate high vibrational levels in the ground electronic state, fluorescence cannot directly 
populate high rotational levels, and no direct paths exist to the levels v" = 3, J" = 25 and 
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v" = 5, J" = 18. H 2 can form in highly excited states (Kokoouline, Greene, & Esry 2001; 
Strasser et al. 2001) via dissociative recombination of fht . Given a dense medium and a 
large ionizing flux, then may be present in large enough quantities to react with H 2 to 
produce Hg in highly excited states. More attempts to detect Hg would be useful. 



7. SUMMARY 

We have modeled a plane-parallel slab of warm H 2 (see Fig. 1), which represents the 
surface of a protoplanetary disk irradiated by UV flux from the star. Our analysis of the 
ultraviolet spectrum of TW Hya obtained with the STIS instrument on HST and with FUSE 
leads to the following conclusions: 

1. The FUV continuum rises shortward of 1700 A, which is indicative of an H 2 dis- 
sociation continuum, although it could also be produced by H 2 fluorescence due to FUV 
continuum pumping or an additional accretion or activity component. The FUV continuum 
is not well explained by simple models of a pure hydrogen slab, which are commonly invoked 
to analyze the excess NUV continuum. 

2. The extinction towards TW Hya is negligible, based on the hydrogen column density 
in our line of sight and assuming an interstellar gas-to-dust ratio. 

3. Self-absorption of Lyman-band transitions involving low excitation energy levels 
weakens the flux in these lines. We model this effect by simulating Lya emission entering a 
plane-parallel atmosphere and pumping the H 2 . 

4. Using the observed H 2 fluxes and our fluorescence models, we reconstruct the Lja 
line profile incident upon the warm molecular layer and compare it to the observed Lja 
profile, for a range of assumed temperatures and column densities. Undetected progressions 
rule out a large region of parameter space in this model. The reconstructed Lja profile is 
similar to the observed profile in the wings, but shows a much narrower absorption feature 
than is observed. This narrow absorption component in the reconstructed profile could be a 
self-reversal or a component of the wind between the source of Lya emission and the warm 
molecular region. 

5. Our models indicate that a molecular layer with a kinetic temperature of 2500lgQQ 
K and a column density of logiV(H 2 ) = 18. (2a error bars) absorbs Lja radiation 
in the surface layer and inner edge of the disk within 2 AU of the central star. The Lya 
pumping leads to a small H 2 dissociation rate and does not cause significant disk dissipation. 
The filling factor of the warm H 2 around TW Hya is rj = 0.25 ± 0.08, although significant 
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uncertainties in the geometry of the fluorescent H 2 weaken our confidence in the large filling 
factor. 

6. The warm H 2 most likely resides in a warm surface layer of the disk. This surface 
layer may be analgous to a classic PDR, although in this case the FUV radiation field is 
dominated by emission lines rather than a continuum. In particular, Lja comprises about 
85% of the FUV radiation field below 2000 A, and it controls the excitation and ionization 
of the disk surface. 

7. Some of the observed H 2 lines are pumped from high rotational levels that cannot 
be excited thermally or by fluorescence. Formation of rotationally excited H 2 by reactions 
with Hjj~ may explain these emission lines. 
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Fig. 1. — A schematic diagram of our model of TW Hya and its disk (not to scale). FUV 
radiation, primarily Lya, is produced at the accretion shock and heats the disk. The warm 
H 2 is located in surface layers of the disk within 2 AU of the star. The hot surface of the 
disk, composed of atoms and ions, probably extends inwards to about 6i?* (Johns-Krull & 
Valenti 2001b), where it is truncated by the stellar magnetosphere. 
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Fig. 2. — The solid line shows the observed spectrum of of TW Hya, after subtracting the 
template spectrum of V819 Tau. The dashed line shows the synthetic spectrum of a pure 
hydrogen slab fit to the observed Balmer jump and the Balmer lines. The pure hydrogen 
continuum accurately approximates the NUV continuum of TW Hya but significantly un- 
derestimates the FUV continuum. The rise of the FUV continuum shortward of 1700 A is 
indicative of an H2 dissociation continuum produced by collisions with energetic electrons, 
although it could also be produced by H2 fluorescence due to FUV pumping or an additional 
accretion or activity component. 
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Fig. 3. — The Lya emission from TW Hya is partially absorbed by H I in the line of sight 
to TW Hya. The interstellar and circumstellar media, centered at about km s -1 , account 
for the red side of the absorption. The blue side of the absporption is produced by H I in 
the wind. The dashed line shows our best fit to the observed Lya profile (solid line) using 
two Gaussian emission profiles to model the intrinsic Lya line before the absorption, while 
the dotted and dashed-dotted lines assume different H I column densities that poorly fit the 
red side of the observed profile. 
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Fig. 4. — The transmission percentage across the observed Lya profile (solid) for column 
densities logiV(H) = 19.75 (dashed) and logiV(H)= 20.5 (dotted), assuming no Doppler 
broadening. This assumption places a strong upper limit of iV(H I< 19.75) towards TW 
Hya. 
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Fig. 5. — The O VI emission line from TW Hya (solid), as observed with the LiFla detector 
on FUSE. We place a limit of iV(H 2 )< 18 in our line of sight towards TW Hya by comparing 
the observed emission with the simulated transmission percentage given the resolution of 
FUSE, through a medium of H 2 with T = 100 K and 7V(H 2 )=17.5 (dashed) and JV(H 2 )=19.0 
(dotted). We would observe many H 2 absorption lines against the O VI emission from TW 
Hya for hotter excitation temperatures of H 2 . 
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Fig. 6. — The NUV continuum (solid) is well fit by the accretion model (dashed) for Ay = 0.0, 
but if Ay = 0.5, the rise in the continuum below 2600 A is not fit well. 
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Fig. 7. — For four different upper levels, we plot the normalized ratio of the observed de- 
excitations in each H 2 line divided by the respective branching ratio (circles). The squares 
show the fluxes prediced based on the optical depth in each lower level for our adopted model 
for logiV^rFj = 18.5 and T = 2500 K (see §5). The dotted and dashed lines show extinction 
curves for Ay = 0.5 and 1.5, respectively. Emission in an optically thin medium, with no 
intervening extinction, would result in ratios of 1.0. The apparent wavelength dependence 
results from the larger optical depth in those levels of H2 which have lower excitation energies. 
The H2 line fluxes from v'—A, J' =18 show no wavelength dependence because the observed 
lines all have lower levels with energies large enough that they are not populated. 
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Fig. 8. — Two plausible geometries for H 2 fluorescence in the surface layers of a protoplane- 
tary disk. The Lja enters the surface layer of H2 at angle 6 with respect to the normal and 
is reprocessed into H 2 emission, which escapes the disk at angle 0. 
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Fig. 9. — 1 —5a confidence contours calculated from Ax„h2 f° r n ^s to the individual H2 line 
fluxes for the thick disk model (see §3). The star indicates the position of the model with 
T = 2500 K and logiV(H 2 ) = 18.5 (see §6). 
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Fig. 10. — Selected regions of the HST/STIS spectrum of TW Hya. The narrow features 
are all H2 lines. The dashed lines show our model fit to the H2 emission. The line at 1333.8 
A is supressed by C II in the wind of TW Hya, which places the H 2 interior to the wind. 
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Fig. 11. — l — 5o" confidence contours calculate from A%^ L (shaded regions) for our fit of 
the reconstructed Lja profile to the red wing of the observed Lya profile, compared with 
1 — 5(7 confidence contours (solid lines) for fits to the observed H 2 fluxes. 
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Fig. 12. — Comparison of the observed (solid line) and calculated (circles) Lya emission 
profiles, for models with (top) log A r (H 2 )=18.5, T = 2500 K, and rj = 0.25 (middle) 
logiV(H 2 ) = 19.0 T = 2000 K, and r] = 0.41; and (bottom) logA^(H 2 ) = 19.0 T = 3000 K, 
and rj = 0.20. The arrows show upper limits to the calculated Lya emission based on un- 
detected H 2 progressions. The top plot is the best fit because of the least scatter in the red 
wing of Lya. The bottom model is unacceptable because several undetected progressions 
(v> = 7, J' = 4 pumped at 1214.4 A and v' = 5, J' = 20 pumped at 1217.7 A) imply up- 
per limits of Lya flux that are significantly lower than the flux reconstructed from nearby 
detected progressions. 




Fig. 13. — Similar contours as in Figure 11. The dots indicate the parameter space that we 
rule out because those parameters predict that lines pumped by 7-4 P(5) should be observed, 
and the diagonal lines indicate the parameter space ruled out because many additional other 
progressions should be observed but are not detected. Filling factor contours (dashed lines) 
of 0.125-4.0 are also shown. A filling factor above 1 is unphysical. The three models with 
reconstructed Lja profiles in Fig. 13 are shown as asterisks. 
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Fig. 14. — The reconstructed Lyct profile at the wavelengths of the pumping transitions 
(diamonds) for T = 2500 K, log iV(H 2 ) = 18.5, and r\ = 0.25 corresponds well to a plausible 
Lja profile (dashed line), constructed using a single Gaussian profile (dotted line) and an 
absorption feature of logiV(H I)= 18.7 blueshifted 90 km s _1 from line center. 
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Fig. 15. — Best-fit models to H 2 fluxes as a function of exit angle of H 2 photons, for 
the thick disk geometry (solid line), thin disk geometry (dashed line), and cloud geometry 
(dotted line). 
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Fig. 16. — Best-fit models to H 2 fluxes as a function of extinction within the H 2 slab, for 
the thick disk geometry (solid line), thin disk geometry (dashed line), and cloud geometry 
(dotted line). The extinction is for a progression in which the average Lya photon is absorbed 
at the midplane of the slab. 
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Fig. 17. — The percent of the slab depth that a Lja photon excites H 2 to two different upper 
levels (absorption depth d a bs in Table 2), at two temperatures. As the column density and 
temperature of the slab increase, the optical depth in the lower levels increase, and a typical 
Lja photon is reprocessed into H 2 at a shallower depth. This transition occurs at different 
column densities and temperatures for each progression, so that the reconstruction of the 
Lja profile is sensitive to these parameters. 
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Table 1: LIST OF OBSERVATIONS 
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1 The FUSE instrument consists of several channels. 
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1 Lower energy level 

2 Radiative decay coefficient (Abgrall et al. 1993) 
3 Absorption depth of Lya photon 
4 Equivalent width of absorption profile within slab 
5 Correction for unseen lines from upper level 

6 Dissociation probability from upper state per electronic excitation, a{b) = a x 10 b . 
7 Dissociation probability from best-fit model. 
8 Average kinetic energy released due to H 2 dissociation. 
9 H 2 dissociation rate for best-fit model (1O _12 M yr _1 ). 
10 Lyo; flux absorbed by lower level (10~ 14 erg cm -2 s _1 ). 
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Table 3: STRONGEST PREDICTED LYMAN H 2 LINES IN THE FUSE BANDPASS 



ID 


^calc 


Emod 


Eobs 


4-0 R(3) 


1053.973 


1.1 


< 3 


4-0 P(5) 


1065.593 


1.6 


< 2 


4-1 R(3) 


1101 889 


2.4 


< 10 


1-0 P(5) 


1109.311 


1.0 


< 6 


3-1 P(l) 


1113.877 


1.4 


< 6 


4-1 P(5) 


1113.949 


1.4 


< 6 


1-1 R(3) 


1148.701 


5.0 


5.5(1.1) 


0-1 R(0) 


1161.693 


1.9 


13.1(b) 


1-1 P(5) 


1161.814 


6.2 


13.1(b) 


1-1 R(6) 


1161 949 


1.8 


13.1(b) 


0-1 R 1 


1162.170 


2.4 


13.1(b) 


4-1 R(12) 


1164.596 


1.1 


< 10 


2-0 P(13) 


1165.834 


1.0 


< 6 


0-1 P(2) 


1166.255 


2.6 


< 6 


0-1 P(3) 


1169.751 


2.5 


< 10 


4-0 R(17) 


1176.325 


1.0 


< 15 


3-1 R(12) 


1179.472 


1.3 


< 10 


1-1 P(8) 


1183.309 


2.7 


< 12 


2-1 R(ll) 


1185.224 


2.8 


< 12 



F obs refers to the observed fluxes and 
Emod to the model fluxes in units of 
10~ 15 erg cm -2 s _1 ,lo" error in (). 
(b) indicates the line is blended with other H 2 lines. 



Table 4: UNDETECTED PROGRESSIONS 



Pump A ca i c A u i E" P Dis 

Em&x A max Kef IDre/ 

A 10~ 7 s _1 eV erg cm -2 s" 1 A A 

1- 1 R(12) 1212.543 4~6 L49 5^5) 3 x 10~ 15 1578 1212.425 1-1 P(ll) 
7-4 P(5) 1214.421 2.8 2.07 0.19 5 x 10" 15 1582 1214.465 3-1 R(15) 
5-3 P(8) 1218.575 6.6 1.89 0.04 5 x 10~ 15 1608 1218.521 2-1 R(14) 
5-0 P(20) 1217.716 5.4 2.39 0.48 2 x HT 15 1223 1217.904 2-1 P(13) 

2- 2 P(8) 1219.154 10.8 1.46 7(-5) 5 x 10~ 15 1579 1219.089 0-2 R(2) 

2- 2 R(9) 1219.101 12.9 1.56 2(-4) 5 x 10" 15 1572 1219.089 0-2 R(2) 
5-3 R(9) 1219.106 5.4 1.99 0.085 5 x 10~ 15 1601 1219.089 0-2 R(2) 

3- 3 R(l) 1215.541 0.46 1.48 2(-4) 5 x 10~ 15 1596 1215.726 1-2 R(6) 
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Table 5: BEST-FIT PARAMETERS FOR ALTERNATE GEOMETRIES 



Geometry" Fit to H 2 line opacities Fit to Lya profile 



Number 


Type 


e 





A V (D) 


iV(H 2 ) 


T 


v b 


V 2 


iV(H 2 ) 


T 


VH2 




1 


Thick 





89 





17.0 


2200 


4.5 


1.58 


18.7 


2400 


0.29 


0.9 


2 


Thick 





77 





17.9 


2600 


0.41 


2.22 


18.5 


2500 


0.25 


2.3 


3 


Thick 





77 


0.5 


17.5 


2700 


2.5 


2.22 


18.8 


2800 


0.23 


9.5 


4 


Thick 





77 


1.0 


17.6 


2300 


16 


3.16 


18.7 


5300 


0.09 


18 


5 


Thick 





77 


1.5 


23.0 


1700 


0.05 


4.21 


19.0 


5300 


0.09 


17 


6 


Thick 





43 





19.1 


2600 


0.12 


4.15 


18.2 


2700 


0.23 


3.0 


7 


Thick 











19.5 


2500 


0.09 


4.98 


18.4 


2600 


0.21 


3.3 


8 


Thin 


80 


89 





16.6 


2600 


1.10 


1.58 


18.3 


2200 


0.31 


1.1 


9 


Thin 


80 


77 





18.7 


2500 


0.15 


3.97 


17.9 


2400 


0.27 


0.8 


10 


Thin 


80 


43 





20.2 


2500 


0.05 


6.94 


18.0 


2400 


0.23 


0.8 


11 


Thin 


80 








20.0 


2500 


0.05 


7.87 


18.0 


2400 


0.22 


0.8 


12 


Thin 


80 





0.5 


19.2 


2500 


0.41 


7.73 


15.3 


3300 


10 5 


29 


13 


Thin 


80 





1.0 


17.9 


1000 


10 13 


6.16 


19.3 


4700 


0.03 


33 


14 


Thin 


80 





1.5 


15.9 


4600 


10 12 


5.86 


19.4 


4900 


0.03 


34 


15 


Cloud 


80 


89 





16.2 


2800 


8.2 


1.96 


18.5 


2400 


0.38 


1.0 


16 


Cloud 


80 


77 





17.3 


2500 


1.2 


2.71 


18.7 


2400 


0.26 


0.5 


17 


Cloud 


80 


43 





17.8 


2500 


0.52 


2.92 


18.7 


2400 


0.25 


1.1 


18 


Cloud 


80 








17.8 


2600 


0.64 


2.80 


18.8 


2300 


0.33 


0.9 


19 


Cloud 


80 





0.5 


18.2 


2000 


12 


2.74 


18.8 


2400 


3.32 


1.4 


20 


Cloud 


80 





1.0 


18.1 


1800 


400 


3.88 


18.8 


2500 


40.4 


1.9 


21 


Cloud 


80 





1.5 


18.1 


1400 


10 4 


6.07 


18.8 


2700 


400 


2.8 



l 9 and are defined in §3.2 and Fig. 8 

'rjH2 calculated from Lya reconstruction for parameters of best fit to H 2 fluxes. 
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Table 6: FLUXES OF STRONG LINES 



ID 


Mab 


Flux 1 




A 


10~ 12 erg cm -2 s 


Total flux 


1170-1700 


65 


Lya 2 


1215.67 


44.74,5 


Lya 3 


1215.67 


80-160 


Continuum 


1170-1700 


6.0 6 


H 2 (obs) 


< 1700 


1.94 10 


H 2 (mod) 


< 1700 


2.20 11 


C IV 


1548 


1.86 6 


He II 


1641 


1.35 


C IV 


1551 


0.98 


C III 


977 


0.44 4 - 8 


C III 


1175 


0.42 


C III 


1175 


0.35 9 


I 


1306 


0.37 


VI 


1032 


0.31 7 


N V 


1238 


0.30 6 


I 


1305 


0.30 4 


I 


1302 


0.28 4 - 5 


VI 


1038 


0.15 7 


C II 


1336 


0.23 4 


N V 


1243 


0.12 


C II 


1335 


0.11 4 - 5 


Si IV 


1393 


0.076 6 


Si IV 


1402 


0.034 6 


III] 


1666 


0.023 



1 These fluxes are more accurate than those in Table 3. 

2 Observed Lya flux. 

3 Flux estimated from reconstructed Lya profile. 
4 Wind absorption attenuates some blue emission. 
5 IS/CS absorption absorbs some emission. 
6 Flux from blended H2 lines subtracted. 
7 Average flux from FUSE LiFlA, LiF2B channels. 
8 Average flux from FUSE SiC2A, SiClB channels. 
9 Average flux from FUSE LiF2A, LiFIB channels. 

10 Observed H 2 flux. 

11 Total H 2 flux, including weak fluxes from model. 



